Distribution of Interacting Ionic Particles in Disordered Media 
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Equilibrium distribution of interacting ionic particles in a charged disordered background is 
studied using the nonlinear Poisson-Boltzmann equation. For an arbitrarily given realization of the 
disorder, an exact solution of the equation is obtained in one dimension using a mapping of the non- 
linear Poisson-Boltzmann equation to a self-consistent Schrodinger equation. The resulting density 
profile shows that the ions are delocalized, despite what the equivalent Schrodinger formulation in 
one dimension would suggest. It is shown that the ions are not distributed so as to locally neutralize 
the background, presumably due to their mutual interactions. 
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Dynamics of interacting charged particles in a hetero- 
geneous Coulombic background, which is a medium with 
a quenched distribution of charge centers, is a complex 
phenomenon involving long equilibration times, a mul- 
titude of time scales and various intermediate regimes, 
and competing mechanisms (Q-^). All these complica- 
tions, which are primarily due to the long-range nature 
of Coulomb interaction, have so far made it extremely dif- 
ficult for even rather sophisticated simulation techniques 
to provide us with a conclusive picture |^|||. 

Perhaps a most preliminary step is to determine 
whether the equilibrium distribution of the mobile ions 
corresponds to a localized or an extended state, in analogy 
to the quantum mechanical (QM) problem of disordered 
electronic systems ||. The answer to this question is 
determined by a delicate balance between the following 
competing effects: (i) attraction of the mobile ions to the 
opposite fixed charges and the tendency towards local 
neutralization, (ii) repulsion between similarly charged 
mobile ions, and (iii) entropy of the mobile ions. (Note 
that for simplicity of the argument we have assumed that 
all the mobile ions are similarly charged whereas all the 
fixed charges are of the opposite type.) While attraction 
tends to localize the mobile ions in the neighborhood of 
the charge centers, repulsion and entropy both favor de- 
localization. One can imagine that determining the out- 
come of this balance in simulations could be a very dif- 
ficult task due to the progressively long escape times for 
the deep Coulomb potential wells. On the other hand, it 
is also not clear if one can trust approximate theoretical 
schemes, which may artificially weaken or strengthen any 
of the above tendencies, in making feasible predictions to 
this end. 

For example, it has been shown that if we neglect 
the mutual repulsion between the mobile ions, the an- 
swer to the above question depends on the dimension- 
ality |||5j@,|8|]. One can show that the disorder-averaged 
density-density correlation function for the mobile ions 

behaves as (c(O)c(x)) ~ exp — (x/C) 2 ~ C 



2 (with £ being the localization length). On general 
grounds one might then expect that the interaction be- 
tween the mobile ions could be effectively taken into 
account by introducing the notion of screening, which 
would change the density-density correlation function to 

<c(0) C (x)> 



exp 



, where g is an ex- 
ponentially decaying function and is the correspond- 
ing Debye screening length that is related to the den- 
sity of the mobile ions. The appearance of the expo- 
nential damping now controls the divergence for d < 2, 
and consequently the mobiles ions will be delocalized 
in any dimension. The simplistic implementation of 
screening using Debye-Hiickel theory, however, is based 
on a linearized approximation of the nonlinear Poisson- 
Boltzmann (PB) theory, which is known to break down in 
the vicinity of the charge centers where the mobile ions 
could presumably be trapped for relatively long times. 
Therefore, it seems that a reliable answer to the local- 
ization question can only be obtained through an exact 
treatment of the full nonlinear PB formulation. This is 
what I attempt to do in this article. 

I consider a one dimensional medium with a quenched 
distribution of charge centers and a finite density of op- 
positely charged mobile ions, and show that one can 
map the corresponding nonlinear PB equation to a self- 
consistent Schrodinger equation. An exact solution of 
the resulting equation is then presented for an arbitrary 
realization of the disorder. The resulting density profile 
shows that the ions are in an extended state. It is shown 
that the mutual interactions of the mobile ions could lead 
to a gross violation of the anticipated local neutrality. 

Consider a set of fixed negatively charged ions de- 
scribed by the given charge density Pf(x), embedded in 
a cloud of positively charged mobile ions with concentra- 
tion c(x, t). The electric potential 0(x,t) is then given by 
the solution of the Poisson equation: 



in d dimcn- 



V 2 0= — (ec + pf), 



(1) 



sions H, and therefore the ions are localized for d < in which e is the dielectric constant of the background, 
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and e is the electron charge. The dynamics of c(x, t) is 
governed by a conservation law of the form dtc+ V- J = 0, 
in which J = —DVc + cv is the current density, where 
v(x, t) is the velocity field corresponding to the mobile 
ions, and D is their diffusion constant. In a mean-field 
approximation, the velocity field of the ions at each point 
is determined by the local value of the electric field E(x, t) 
as v = /iE = — /iV0, where p is the electric mobility of 
the ions. Using the self-consistency relation, the conti- 
nuity equation can be written as 



d t c = DV 2 c + pV ■ (cV(/>). 



(2) 



The above two equations describe the nonlinear dynam- 
ics of mobile ions in the presence of a given distribution 
of fixed ions At equilibrium we have dtc = 0, 

and Eq.(^) yields c ~ exp(— fi(f)/D), where equilibra- 
tion at a temperature T implies an Einstein relation 
D = (n/e)kBT. Inserting the Boltzmann form for c back 
into the Poisson equation then yields the PB equation. 

Let us now restrict ourselves to the one dimensional 
case where the equations prove to be more tractable. 
Eliminating the concentration field c from the two equa- 
tions, one obtains an equation for the electrostatic po- 
tential, which can be integrated twice to yield 



d t <j>{x, t) - ft^-oo, t) = Dd 2 x <j) + |(a^) 2 



(3) 
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where the boundary condition _E(±oo,i) = has been 
implemented. The above equation can be called a dy- 
namical PB equation Jl2| . 

We focus, for simplicity, on the specific example of a 
series of TV + 1 negatively charged permeable plates with 
charge (number) density a and surface area A, which 
are placed at random locations x n on the a;-axis, i.e., 

p f {x) = ~eaJ2n=o 6 ( x ~ X n)- 

Using a so-called Cole-Hopf transformation <f){x,t) — 
cf)(—oo,t) — (2D//i) In W(x,t), Eq.(||) can be written as 
a Schrodinger-like equation for W(x, t): 



d t W(x,t) = Dd 2 x W(x,t) 

N 



(4) 



2D 



[e n (t)9(x - x n ) - S(x - x n )} W(x, t), 



n=0 



where 



.(*) 



d x W(x+,t)+d x W(x-,t) 
W(x n ,t) 



(5) 



and A = eD/(ireap) = eksT/ (ire 2 a) is a characteristic 
length scale corresponding to the relative strength of the 
electrostatic and thermal interactions. Note that Eq.(|J) 
is still nonlinear due to the self-consistent coupling of the 
so-called asymmetry fields e n (t) = pE(x n ,t)/D. 

The value of the asymmetry field at each point can be 
calculated from Gauss' law as 



r*\ 2 ( 2Q A 



(6) 



in which Q n is the net amount of (uncompensated) charge 
on the left hand side of point x n - (Overall neutrality then 
implies that — Q n should be the net amount of charge on 
the right hand side of point x n .) In other words, s n (t) 
measures the relative deviation from local neutrality as 
seen at x n - 

The dynamics described by Eq.(|J) is difficult due to 
the time dependence of the e n (t) fields. We thus focus 
on the long time limit of the dynamics, when the ions 
are expected to equilibrate in the disordered medium. In 
this limit, the asymmetry fields are time independent, 
and the Schrodinger equation reads 



A 



d 2 x W{x) + V(x)W(x) = 0, 



where the effective potential is given as 



N 



V(x) = 2J [-£n®(x - X n ) + S(x - X n )} . 



(7) 



(8) 



Note the peculiar feature that the delta function peaks 
in the potential, which correspond to the attractive sites 
for the mobile ions, appear as effective repulsive potential 
wells in the equivalent QM description! 

To complete the QM formulation, we need a nor- 
malization for W(x). Using the equilibrium relation 
c ~ exp(— /j,<p/D), and the definition of W, we can re- 
late c to W as c(x) ~ 1/W 2 (x). Conservation of the 
overall number of mobile ions can thus be written as 
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W{x)_ 



(9) 



where L is the system size. 

We can now solve the nonlinear self-consistent QM 
problem described by Eqs.(||), (0), and ([))), because the 
potential is a constant in each interval. The solution of 
Eq.(^) is a series of plane waves: 



W(x) 



where 



A Q + B x, 



Aj e ik i( x - x i) 



for x < Xq , 



+Bj e ik i( x x i\ for xj^i < x < 
> A N+ i + B N+ i x, for x N < x, 



(10) 



k 3 = 



j'-i 



n=0 



1/2 



(11) 



for j = 1, • • ■ , N. We should then match the neighbor- 
ing solutions across the boundary points based on the 
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following three requirements: (i) W(x) should be con- 
tinuous, (ii) the jump in d x W/W should be set to 2/A 
as implied by the delta function peaks in V(x), and (hi) 
the asymmetry fields should satisfy the self-consistency 
relation of Eq.(||). These criteria translate into a series 
of self-consistency equations for the wave vectors: 
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(12) 



for j = 1, ■ • • ,N, where Aj = Xj —Xj-i and fco = fcjv+i = 
0. For any given set of AjS, one should solve these cou- 
pled nonlinear equations to determine the wave vectors. 
The matching requirements together with the normaliza- 
tion condition of Eq. (||) , can then give us the full solution, 
and thus the density profile of the ions. 

In order to determine whether the ions are in an ex- 
tended state, or localized near the adsorbing sites, we 
need to find out whether the fcjS that satisfy Eq.(|l2|) are 
real or imaginary. Instead of trying to solve for fcjS ex- 
actly, I shall first resort to approximate schemes. Physi- 
cally we expect to have a tendency towards localization 
at low temperatures or for strong disorder, and nearly 
free ions with uniform distribution at high temperatures 
or for weak disorder. 

Let us first focus on the case of strong disorder where 
we have A < A 3 -, and try to solve Eq. ([12]) perturbatively. 
I find 
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for the first few terms, which is real. This proves, strictly 
speaking, that if the ions are localized, the localization 
length should be greater than A 5 /A 4 . However, it looks 
rather unlikely that higher order corrections start to de- 
velop imaginary components, which implies that the ions 
should be delocalized even in the strong disorder limit. 

To check this, I have solved Eq. (|l2|) numerically for 
various choices of {Aj}, and the results for the values of 
{kj} were always real. For example, for the case of 11 
charged centers with Ai/ A = 0.1, A 2 /A = 0.5, A 3 /A = 2, 
A 4 /A = 0.8, A 5 /A = 0.01, A 6 /A = 7, A 7 /A = 3.2, 
A 8 /A = 0.3, A 9 /A = 1.6, and Am/A = 0.0002, I find 
fciA = 33.84, k 2 X = 7.33, fc 3 A = 1.42, fc 4 A = 4.95, 
fc 5 A = 316.69, k 6 X = 0.399, k 7 X = 0.757, k$X = 12.599, 
kgX — 1.933, and fcioA = 15710, which are all real. This 
confirms that the ions are in an extended state. 

We can also calculate the asymmetry fields from the 
above results as 



to the leading order. This shows that unlike what one 
might expect, the distribution of ions is not predomi- 
nantly determined by the criterion of local neutralization 
of the background, which would yield Ej =0. 

Now let us address the peculiar features of the for- 
mulation, and try to understand them better: (i) Why 
do the delta function peaks, corresponding to adsorbing 
sites in the real problem, appear as repulsive in the effec- 
tive potential V{x)l (ii) While we expect a Schrodinger 
equation with random potential in one dimension to lead 
to localization H, why do we find that the ions are de- 
localized? To answer these questions, we should note 
the unusual normalization condition of Eq.(|^) above. A 
true wave function for this problem should be defined as 
9(x) — 1/W(x) according to Eq.(j|), leading to a corre- 
sponding Schrodinger equation as 



A 



d x V(x) 



V(x) 



= 0, (15) 



which is highly nonlinear. 

While the sign of the effective potential is now reversed 
corresponding to attractive delta function peaks, an ad- 
ditional effective repulsive self-consistent potential has 
also appeared. This potential, which is proportional to 
the the electrostatic pressure ~ (Vtfi) 2 , disfavors pile- 
up of the ions in the vicinity of the attractive sites and 
promotes derealization. This extra self-consistent term, 
which is a direct consequence of the repulsive interac- 
tions between the mobile ions, explains the discrepancy 
between the predictions of this theory and conventional 
Schrodinger equation in a random potential. 

In conclusion, I have presented an exact analysis of 
the problem of interacting mobile ions in a random back- 
ground of fixed ions in one dimension. The density pro- 
file of the ions at equilibrium is shown to be extended 
(as opposed to localized) as a consequence of the interac- 
tions between the mobile ions. Since the one dimensional 
problem in the strong disorder limit that has the best 
chance of leading to localization is proven to result in 
an extended state, it is plausible to assert that the same 
result should hold true in general, i.e., in higher dimen- 
sionalities and for weaker disorder. Although I have only 
considered positively charged mobile ions in a negative 
random background, the formulation can be readily gen- 
eralized to contain both species of mobile ions and fixed 
ions. 
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